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Abstract 

We discuss the intermediate wave-packet formalism for analytically quantifying the energy de- 
pendence of the two-flavor conversion formula that is usually considered for analyzing neutrino 
oscillations and adjusting the focusing horn, target position and/or detector location of some fla- 
vor conversion experiments. Following a sequence of analytical approximations where we consider 
the second-order corrections in a power series expansion of the energy, we point out a residual time- 
dependent phase which, in addition to some well known wave-packet effects, can subtly modify the 
oscillation parameters and limits. In the present precision era of neutrino oscillation experiments 
where higher precision measurements are required, we quantify some small corrections in neutrino 
flavor conversion formulae which lead to a modified energy-dependence for z/„ «-> u e oscillations. 
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Although neutrino physics has been fueled by the recent growth in quality and quantity of 
experimental data Q, Q, E , 3, 3,3], there are still open questions on the theoretical front 
3, Q, 14, 15] which, in some cases, concern with obtaining more refined parameters from 



flavor conversion formulae 0, 0, Q, Q, |20| . In the current experimental scenario, we can 
notice that the KamLAND experiment will significantly reduce the allowed region for Am 2 2 
and sin (29 12 ) parameters, where the second-order wave-packet corrections can appear as an 
additional ingredient for accurately applying the phenomenological analysis 2l|. In parallel, 



the next major goal for the reactor neutrino program will be to attempt a measurement of 
sin 2 (29 13 ), i.e. while the determination of the mixing parameters appearing in the solar 

HQ 

and atmospheric [3,0, [5] neutrino oscillations has already entered the precision era, the next 
question which can be approached experimentally is that one of e3 mixing. Some important 
experiments which will search for more precise measurements of 9 13 are Double Chooz [k]] 
with designed sensitivity of 0.03, and the Daya Bay Reactor Neutrino Experiment [ll| with 
designed sensitivity of 0.01. A consequence of a non-zero U e3 matrix element will be a small 
appearance of v e in a bean of v^. Assuming the scenario where Am 2 2 -C Am 2 3 , which is 
suggested by experimental data, and for E v ~ L Am 2 3 /(27r), ignoring matter effects, we find 



Piffj, —> v e ; L, E v ) ~ sin 2 (26*13) sin 2 (#23) sin 2 



Am 2 23 L 



±E V 



(1) 



This expression illustrates that 6 13 manifests itself in the amplitude of an oscillation between 
the second and third families. To improve the experimental limits on 8 13 , one needs both 
good statistics and low background data. At the same time, all kind of fine-tuning correction 
should also deserve a quantitative analysis. In particular, it can be shown that reactor 
experiments have the potential to determine 9 13 without ambiguity from CP violation or 
matter effects (by assuming the necessary statistical precision which requires large reactor 
power and large detector size). With reasonable systematic errors (< 1%) the sensitivity is 
supposed to reach sin 2 (26\ 3 ) « 0.01 — 0.02 [l^] and an accurate method of analysis, maybe 
in the wave-packet framework, can be required. 

The most simplified theoretical formulation used for describing the flavor conversion 



process involves the intermediate wave-packet treatment 



22 



231 ] which eliminates the most 



controversial points rising up with the standard plane-wave formalism 



22 



24 



25|[34|. It is 



convenient to observe that the intermediate wave-packet procedure leads to flavor conversion 
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expressions that, after some suitable parameter adjustments, are mathematically equivalent 



to those ones obtained with the average energy treatment usually applied to plane waves [26]. 
From the point of view of a first quantized theory and in the context of vacuum oscillations, 
our main purpose is to compare the standard quantum oscillation plane wave treatment with 
an analytical study in the wave-packet framework by re-obtaining the energy dependence of 
the oscillation probability formula in a particular phenomeno logical context. In this sense, 
by analytically quantifying the dependence of the neutrino oscillation parameters on the 
product between the wave-packet width a and the average energy e of detection, we shall 
obtain the their range of deviation from the values obtained with the plane-wave approach. 
Therefore, we suggest an improvement on bounds in adjusting the focusing horn, target 
position and/or detector location for some flavor conversion experiments. 

In neutrino oscillation experiments, the distance of the detector from the source L, the 
neutrino average energy E v = e, and the appearance (disappearance) probability (P) are the 
experimental input parameters which lead to the output mixing angle and mass-difference 
parameters. For discussing how the procedure for obtaining this parameters can be modi- 
fied/improved, we are effectively interested in quantifying the energy dependence of oscilla- 
tion probabilities when the intermediate wave-packet treatment is taken into account. 

The first step of our study concerns the analytical derivation of a flavor oscillation for- 
mula where a gaussian momentum distribution and a power series expansion of the energy 
up to the second-order terms are utilized for obtaining analytically integrable expressions 
which result in the flavor conversion probabilities. We also state that the main aspects of the 
oscillation phenomena can be understood by studying the two-flavor problem. In addition, 
substantial mathematical simplifications result from the assumption that the space depen- 
dence of wave functions is one-dimensional (z-axis). With such simplifying hypotheses, the 
time evolution of flavor wave-packets can be described by 

$>(z,t) = 4> 1 (z,t) cos 9 V! + 4> 2 (z , t) sin 6 u 2 

= [<j>x{z,t) cos 2 9 + 4> 2 {z,t) sin 2 9) v a + 

[4> 2 (z,t) — cp 1 (z,t)] cos 9 sin 9 up 

= 4> a (z,t;9) v a + <j)p(z, t; 9) up, (2) 

where u a and vp are flavor-eigenstates and v r and i/ 2 are mass-eigenstates. The probability 
that neutrinos originally created cis db i^ct 

flavor-eigenstate with average energy e oscillate 



into a up flavor-eigenstate after a time t is given by the v@ coefficient 

sin 2 (20) 



dz \<j)p\' 



{ 1 - iNT(t) } 



where iNT(t) represents the mass-eigenstate interference term given by 



iNT(t) = Re 



dz 4>\(z,t) 4> 2 {z,t) 



(3) 



(4) 



Since the time evolution of each mass-eigenstate wave-packet is given by [2J 



<t>i{z,t) 



<P[Pz-Pi)exp 



: t 



l t + ip z z 



where Ep} = (p 2 z + m 2 ) 2 , i = 1, 2 and tp(p z -Pi) = (27ra 2 ) 4 exp 
the interference term iNT(i) by evaluating the following integral 

dp 



{Pz-Vi) 2 a 2 



(5) 

we can calculate 



/ 



2vr 



<p(jPz-Pi)<p(pz-Pi) exp [-i AE Pz t] 



exp 



-(aAp) 2 



dp z 



2tt 



<P 2 (Pz-Po) exp [-i AE Pz t], 



(6) 



where we have changed the ^-integration into a p z -integration and introduced the quantities 
Ap = Pi-p 2 , p = \{Pi + Pz) an d AE Pz = E p ] -E p ] . The oscillation term is delimited by 
the exponential function of aAp at any instant of time. Under this condition, we could 
never observe a pure flavor-eigenstate. Besides, oscillations are considerably suppressed if 
aAp > 1. A necessary condition to observe oscillations is that aAp 1. This constraint 
can also be expressed by 5p 2> Ap where Sp is the momentum uncertainty of the particle. 
The overlap between the momentum distributions is indeed relevant only for 5p ^> Ap. 
Strictly speaking, we are assuming that the oscillation length (7r . 4e 2 ) is sufficiently larger 
than the wave-packet width, which simply says that the wave-packet must_not be as wide as 
the oscillation length, otherwise the oscillations are washed out 
to the Eq. (jBJ), without loss of generality, we can assume 
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31 



321 ] . Turning back 



iNT(i) =Relj P 2 (P*-Po) exp [-i &E Pz t] | 



(7) 



In the literature, this equation is often obtained by assuming two mass-eigenstate wave- 
packets described by the same momentum distribution centered around the average momen- 
tum p = p . This simplifying hypothesis also guarantees the instantaneous creation of a pure 

In fact, we get (f> a (z, 0, 6) = (f>i(z, 0) = 4> 2 (z, 0) from Eq. (J2J) 



flavor eigenstate u a at t = 



and (j)i3(z,Q,6) = 0. In order to obtain an expression for (pi(z,t) by analytically evaluating 



the integral in Eq. (j5J) we firstly rewrite the energy Ep] as Ep* = E 4 [1 + a t (a t + 2v i )] 5 , 
where E 4 = (m\ + p 2 )^ , v 4 = |r and = p " ~ Po . The integral in Eq. (J5]) can be analytically 
evaluated only if we consider terms up to order of in a power series expansion conveniently 
truncated as 



HO 



E, 



07 



1 + + y (1 



■m: 



£ i +p oO -, + ^- O ? + 0((7?) 



(8) 



The zero-order term in the above expansion gives the standard plane-wave oscillation 
phase. The first-order term p <7 4 is responsible for the slippage between the mass-eigenstate 
wave-packets due to their different group velocities. It represents a linear correction to the 

2 

standard oscillation phase. Finally, the second-order term ^hof, which is a (quadratic) 
secondary correction, gives the well known spreading effects in the time propagation of 
the wave-packet. Moreover, it leads to the appearance of an additional time-dependent 
phase in the final expression for the oscillation probability. In case of gaussian momentum 
distributions, all these terms can be analytically quantified. By evaluating the integral (H) 
with the approximation (JSJ), and after performing some mathematical manipulations 271. |28|. 
we can express the interference term as 

iNT(t) = DMP(i) X OSC(t), 



which was factorized into a decoherence damping term given by 

(Avt) 2 



2a 2 [l + SP 2 (t)] 



DMP(i) = [1 + Sp 2 (t)]"3exp 
and a time-oscillating flavor conversion term given by 

Osc(t) = Re{exp[-iAEt-iO(t)}} 



(9) 



(10) 



cos [AEt + e(t% 



where 



and 



©(*) 



, . t . ( m 2 



\ arctan [Sp(i)] 



a 2 p 



,2 2 

Po 



Sp 3 (t) 



2 p' [i+$r\t) 



(11) 



(12) 



(13) 



with p = 1 — [3 + (^r) 2 ] fr, A£7 = E 1 — E 2 and e = \ / E 1 E 2 . The time-dependent quantities 
Sp(£) and G(t) carry the second-order corrections and, consequently, the spreading effect to 



the oscillation probability formula [29l . l3Qf | . If AE <C e, the parameter p is limited by 

2 

the interval [1, —2] and it assumes the zero value when ^ « |. The slippage between the 
mass-eigenstate wave-packets is quantified by the vanishing behavior of the damping term 
DMP(t). The NR limit is obtained by setting p 2 = 1 and p = in Eq. ( TTOT) . In the same 
way, the UR limit is obtained by setting p 2 = 4 and p = e. In fact, the minimal influence 

2 

due to second-order corrections occurs when f£ ~ 3 (p ~ 0). Returning to the exponential 
term of Eq. ([TO]) , we observe that the oscillation amplitude is more relevant when Avt <C a. 
It characterizes the minimal slippage between the mass-eigenstate wave-packets which occur 
when the complete spatial intersection between themselves starts to diminish during the 
time evolution. 

The oscillating component Osc(t) of the interference term lNT(t) differs from the standard 
oscillating term cos[AEt] by the presence of the residual phase ®(t), which is essentially 
a second-order correction {29], 3(|. Superposing the effects of Dmp(£) and the oscillating 
character Osc(t), we immediately obtain the flavor oscillation probability in its explicit 
form 



27 



P{u a - v#t) « {l - [1 + Sp 2 (t)H 



x exp 



(Avf) 2 
2a 2 [l+Sp 2 (t) 



cos[AEt + e(t)] }, (14) 



from which we notice that the larger is the value of ae, the smaller are the wave-packet 
effects. 

To perform some phenomenological analysis involving the 13 mixing angle, we replace 
sin 2 (20) by sin 2 (20 13 ) sin 2 (0 23 ) and we set AE = AE 23 in the Eq. (|14p in order to realistically 
characterize the — > v e conversion which emerges in a three flavor scenario. We establish 
the experimental input parameters as being the distance L, the neutrino energy distribution 
e and the appearance (disappearance) probability (P). At this point, it is instructive to 
redefine some parameters which shall carry the main physical information in the oscillation 
formula. Firstly, we set the oscillation length scale L which is related to an energy scale 
E by the expression L = 27r A ^ . Both parameters, L and E , correspond to referential 
scales that can be calibrated in agreement with the experimental configuration and the data 
analyzed. From a practical point of view, the criteria for the choice of this parameters is not 
so arbitrary. In a phenomenological analysis, the choice of the parameter E can be done 
in correspondence with the peak of an energy distribution (e) of a certain type of neutrino 
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flux for which the experimentally obtained energy distribution is typically determined by 
the neutrino production processes. As we shall observe in the analysis which follows the 
calculations, in order to quantify the corrections due to the wave packet approximation, the 
reference value of E can also be set equal to an averaged value (e) = E where, depending 
on the width of the energy distribution of the neutrino flux, the necessity of an additional 
energy integration over e (averaged out integration) is discarded. 

We also introduce the auxiliary definitions 5 = as and v — ^ which respectively param- 
eterize the wave-packet character and the propagation regime. With the previous definition 
of £, we introduce the dimensionless variables, 



e 



and I 



L 



(15) 



which will be useful in the the subsequent analysis, since it allows us to extend the validity of 
the interpreted results to any set of parameters L and E . In real experiments the neutrino 
energy, e, and sometimes the detection position, vt ~ L, can have some spread around 
and/or deviation from respectively E and L due to various effects, but in a subset of this 
experiments there is a well-defined value of (L/e) ~ L /E (or x/l ~ 1 in the plane- wave 
limit as we shall see in the following) around which the events distribute. Following the 
same approach that we have adopted while we were analyzing the parameter p in Eq. (fl~2l) . 
if AE <C e, which is perfectly acceptable from the experimental point of view, we can write 
e = y/E 1 E 2 ~ + E 2 ) so that an effective plane- wave flavor conversion formula can be 
obtained from Eq. ([1]) as 



| 1 — COS 


"vrf 


1 








X 





(16) 



Analogously, we can observe by means of the Eqs (191 113]) that the wave-packet flavor con- 
version formula with second-order corrections (fT^|) exhibits a similar implicit dependence on 
time. The Eq. ( TH|) can thus be rewritten as a function of the above parameter x ffTBT) in 
terms of 

Dmp(x) = [1 + Sp 2 (x) 



and 



Osc(x) 



cxp 



cos 



7Tt 



25xvJ 1 + Sp 2 (x) 



+ 6(x) 



(17) 



(18) 
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where 



and 



e(x) 



Sp(x) 



■ arctan [Sp(x) 



Til 

x 5 2 
5 2 v 2 



Sp 3 (x) 



(19) 



(20) 



2 L ~" v /J 2p 2 [1 + Sp 2 (x) 

with p « 1 — 3?A We can observe that the parameter 5 = ae carries the relevant information 
about the wave-packet width a and the average energy e. If it was sufficiently large (5> 1) 
so that we could ignore the second-order corrections in Eq. (jSJ), the probability with the 
leading terms could be read as 



(P) WP1 = sin 2 (20i 3 ) sin 2 

~7lT 



1 — COS 



cxp 



7rZ 



(21) 



which, in the particular case of an ultra- relativistic propagation (v = 1), can be used as 
a reference for comparison with experimental data 2l|. By the way, despite the relevant 
dependence on the propagation regime (v), once we are interested in some realistic physical 
situations, the following analysis will be limited to the ultra-relativistic propagation regime 
corresponding to the effective neutrino energy of the current flavor oscillation experiments. 

As previously mentioned, the shape of the oscillation probability curve as a function 
of the energy (x) for the above approximations is, indeed, different from that one of the 
standard plane-wave treatment, as we can observe in the Fig. ([1]) where we have plotted 
the fixed-distance probabilities P (u^ — > u e ) normalized by sin 2 (20 13 ) as a function of the 
dimensionless energy x — J- for four different values of 5 = ae. In the first plot we 
illustrate the wave-packet approximation with 1 st order corrections parameterized by the 
Eq. ( 12T1) and in the second plot we illustrate to the wave-packet approximation with 2 nd 
order corrections parameterized by the Eq. (11411 where the dependence on x is implicit. The 
quoted approximations can be compared with the plane- wave approximation (dotted line). 
In order to keep clear the meaning of the deviation of the wave-packet approximation from 
the plane-wave approximation, in spite of the dependence on the energy of the parameter 
5 = ae, we are constrained to set constant values to it for each curve which expresses the 
probability dependence on the energy. Alternatively, we could set S = aE and 5 = x 5 in 
order to re-plot the oscillation probability dependence on x, which is, however, completely 
unrealistic under the point of view of the approximation accurateness. The correction on 
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FIG. 1: Fixed-distance probabilities Piy^ — ► u e ) at Z = = 1, normalized by sin 2 (2# 13 ), with 
^23 = tt/4, as a function of the dimensionless energy x = -J^. We have assumed that the wave-packet 
approximation is fixed by 5 = ae = 2, 5, 10, 20 for any arbitrary value of x{e). For sufficiently large 
values of S, for instance when 5 = ae = 20, we recover the plane wave result which, at first glance 
(visually), coincides with the I s * (first plot) 2 nd (second plot) order results respectively given by 
Eqs. (|2ip and (|14p . Just for completeness, the same map can be reproduced when we set x — > x/l 
for unconstrained Z, i.e. instead of assuming Z = 1 which sets the values of L and E separately, 
we can choose to fix the ratio {L /E ) eliminating one degree of freedom. It allows us to extend 
the validity of the information that we can extract from the figure to a larger set of parameters L 
and E a which characterizes an experimental apparatus. 

the first maximum of probability that allows us to adjust the focusing horn, target position 
and/or detector location for some flavor conversion experiments is represented in the Fig. [2] 
where the maximal values of x were numerically obtained as a function of ae. Considering 
the energy dependence represented in the Fig. ([!]), it is advantageous to introduce a third axis 
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representing the dependence on the parameter ae in order to illustrate the complete/effective 
oscillation behavior. The Fig. ([3]) allows us to qualitatively identify the influence of the wave- 
packet corrections brought up by ae. Quite generally, the complete analysis of the oscillating 
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as 



FIG. 2: The first maximum x Max of the probability expression Piy^ — > v e ) as a function of the 
wave-packet correction parameter ae varying from 1 to 20. For plane waves the maximum occurs 



at Xi\ 



1. 



character coupled to the loss of coherence between the mass-eigenstate wave-packets, which 
suppresses the flavor oscillation amplitude, depends on the experimental features such as 
the size of the source, which allows estimating the wave-packet width (a), the neutrino 
energy distribution (e), and the detector resolution (L ). Once they produce an effect 
competing with that of the finite size of the wave-packet, the neutrino energy measurements 
cannot be performed very precisely. If we set the energy uncertainty represented by SE, the 
Heisenberg uncertainty relation states that SE a ~ 1 and, consequently, our approximation 
hypothesis leads to ^ ~ ^ <C 1. Realistically speaking, a typical neutrino-oscillation 
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The Oscillation Probability — > 

sin 2 (2 13 ) 




0.2 0.4 0.6 0.8 1 1.2 1.4 



x 

FIG. 3: The contour plot corresponding to the 3-dimensional representation of the dependence of 
fixed-distance probabilities P{vn — > v e ) (at I = = 1) on the parameter ae which characterizes 
the dependence on the wave-packet approximation. This plot illustrates the magnitude of the 
oscillations in a map of x(e) by ae. Darker corresponds to small probabilities. 

experiment searches for flavor conversions by means of an apparatus which, apart from the 
details inherent to the physical process, provides an indirect measurement of the neutrino 
energy (in each event) accompanied by an experimental error AE exp due to the "detector 
resolution". In case of AE £ xp < ~ ~ -L the effective role of the second-order corrections 

E E aE' 

illustrated in this analysis can be relevant since, as we have anticipated, the necessity of 
an additional energy integration over the energy distribution e (averaged out integration) is 
discarded. On the contrary, AE exp > 5E demands for an average energy integration where 
the decoherence effect through imperfect neutrino energy measurements by far dominates. 
In this sense, the current experimental values/measurements set some limitations on the 
applicability of our analysis which, at this point, is restricted to the 7 Be and pep lines for 
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solar neutrinos ( E £ xp C 1), certainly to some (next generation) reactor experiment where 
the designed sensitivity is of the order of (ae)' 1 , and eventually to supernova neutrinos 35]. 

Generically speaking, although the higher energy neutrinos are more accessible exper- 
imentally, the corrections to the wave packet formalism can be physically relevant for 
p — p neutrinos with energy distributed around the values of e ~ 10 — 100 keV. Follow- 



ing the standard procedure 



26] (which is not free of controversial criticisms) for calcu- 



lating the wave packet width a of the neutrino flux for p — p solar reactions, we obtain 
a = 10~ 10 — 10~ 8 m = 0.5 — 50 (keV)~ l . Such an interval sets a very particular range for the 
a parameter comprised by the interval 5 ■ 10~ 5 — 0.2 for e ~ 0.01 — 0.4 MeV which introduces 
the possibility for wave packet second-order corrections establish some not ignoble modifi- 
cations for the p — p neutrino oscillation parameter limits. In a supernova, the size of the 
wave packet is determined by the region of production (plasma), due to a process known as 
pressure broadening, which depends on the temperature, the plasma density, the mean free 
path of the neutrino producing particle and its mean termal velocity 26j. Neutrinos from 



supernova core with 100 MeV energy have a wave packet size varying from ~ 5 • 10~ 16 m to 
~ 10~ 14 m which leads to a wave packet parameter ae comprised by the interval 0.25 — 5 
for which the second-order corrections can be indeed relevant. In fact, once we have precise 
values for the input parameters (P), L and e, we could determine the effectiveness of the 
first /second-order corrections in determining Am 2 for any class of neutrino oscillation ex- 
periment. For instance, the flux of atmospheric neutrinos produced by collisions of cosmic 
rays (which are mostly protons) with the upper atmosphere is measured by experiments 
prepared for observing v e «-> and v e <-> conversions. The neutrino energies range about 
from 0.1 GeV to 100 GeV which constrains the relevance of WP effects to an wave packet 
width a ~ 10~ 12 m. The majority of the old generation of the reactor neutrino experiments, 
cover a large variety of neutrino flavor conversions where the neutrino energy flux times the 
corresponding wave packet width a makes the wave packet second-order corrections, at first 
glance, not so relevant (ae » 1 tends to the plane wave limit). 

As an additional remark, it is pertinent to emphasize that there is no accurate way to 
experimentally measure or phenomenologically compute the wave-packet width of a certain 
type of neutrino flux, for which we have only crude estimations. Consequently, apart from 



the obvious criticisms to the plane-wave approach 



221 ]. we cannot arbitrarily assume that 



the modifications introduced by the wave-packet treatment (in particular, with second- 
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order corrections) are irrelevant in the analysis of any generic class of neutrino experimental 
data. Maybe, in a very particular scenario, the above analysis can be applied in designs 
of some experiment dedicated to the 9 13 mixing angle measurement. Finally, from the 



phenomenological point of view, the general arguments presented in [21| continue to be 
valid, i.e. the above discussion has so far been limited to vacuum oscillations. In conclusion, 
the characterization of the wave-packet (a) (implicitly described by a) accompanied by the 
precise determination of the neutrino energy distribution (e) should be considered when the 
accuracy in obtaining the neutrino oscillation parameters or their limits is the subject of the 
phenomenological analysis. 
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